arXiv:1506.07554v2 [q-fin.CP] 1 Jul 2015 


Double-jump stochastic volatility model for VIX: 

evidence from VVIX 


Xin Zang* 
Peking University 

Jun Ni^ 

Penn State University 

Jing-Zhi Huang^ 
Penn State University 

Lan Wu§ 

Peking University 


This version: October 31, 2016 


* School of Mathematical Sciences, Peking University, Beijing 100871, PR China, USA; E-mail: 
xzang@pku.edu.cn. 

^Department of Mathematics, Pennsylvania State University, University Park, PA 16802, USA; E-mail: 
jznl32@psu.edu. 

^Department of Finance, The Smeal College of Business, Pennsylvania State University, University Park, 
PA 16802, USA; E-mail: jxh56@psu.edu. 

§Key Laboratory of Mathematical Economics and Quantitative Finance, Peking University, Beijing 
100871, PR China; E-mail: lwu@pku.edu.cn. 


1 



Double-jump stochastic volatility model for VIX: evidence from 

VVIX 


Abstract 

The paper studies the continuous-time dynamics of VIX with stochastic volatility and 
jumps in VIX and volatility. Built on the general parametric affine model with stochastic 
volatility and jump in logarithm of VIX, we derive a linear relation between the stochastic 
volatility factor and VVIX index. We detect the existence of co-jump of VIX and VVIX and 
put forward a double-jump stochastic volatility model for VIX through its joint property with 
VVIX. With VVIX index as a proxy for the stochastic volatility, we use MCMC method to 
estimate the dynamics of VIX. Comparing nested models on VIX, we show the jump in 
VIX and the volatility factor is statistically signihcant. The jump intensity is also state- 
dependent. We analyze the impact of jump factor on the VIX dynamics. 

Keywords: Volatility indices. Volatility proxy. Co-jump, Monte Carlo Markov chain, Bayesian 
analysis 



1 Introduction 


Modelling VIX index and its derivatives has been a hot topic among researchers. As a 
measure for market’s expectation of 30-day implied volatility of S&P500 index, VIX provides 
rich information for the prediction of market’s future trend. It can be seen as a compression 
of information involved in S&P500 options. Usually, the VIX and S&P500 index have an 
empirical negative correlation relationship, so VIX index is often referred to as the fear index 


or the fear gauge. For more about VIX, see e.g. Carr and Wu (2005). 

Out of its importance, much attention has been focused on modelling the dynamics of 
VIX directly. Earlier work tries using geometric Brownian motion, square root diffusion or 
log-normal Ornstein-Uhlenbeck (OU) diffusion to model VIX. Jumps in VIX are also added 
by some authors. Recently, a novel parameterized stochastic volatility model for VIX is put 


forward by Mencia and Sentana (2013) and Kaeck and Alexander (2013). They specify a 
new process to model the volatility of VIX which may be correlated with VIX and show its 
empirical advantage over the other traditional models. Both of them also point out that the 
appearance of this stochastic volatility reduces the impact of jump on VIX. However, the 


model specihcation of this novel stochastic volatility diverges between them. Mencia and 


Sentana (2013) adopt a pure-jump OU process out of the analytical treatability while Kaeck 


and Alexander (2013) characterize the volatility as a square root diffusion which takes the 


correlation of the volatility and VIX into account. How to specify and estimate this volatility 
factor and further interpret the dynamics of VIX better is still an open problem. 

In 2012 OBOE introduces a new volatility index named VVIX into market. Like the role 


of VIX, VVIX measures the 30-day implied volatility of VIX index. Huang and Shaliastovich 


(2014) construct the realized volatility of VIX index (i.e. realized volatility of volatility) 
and show that the VIX index itself is not a good predictor for its realized volatility while 
VVIX serves as a better candidate. We thus can infer from their empirical conclusion that 
the VVIX index may provide some extra information about the volatility of VIX beyond the 
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VIX itself. 


In this paper we mainly study the dynamics of VIX and especially concentrate on mod¬ 
elling its stochastic volatility under the physical measure via additional information provided 
from VVIX index. Based on joint behavior of VIX and VVIX we put forward a double-jump 
stochastic volatility model for logVIX and its volatility. We will show that under a gen¬ 
eral affine assumption for the logarithm of VIX and its stochastic volatility, the stochastic 
volatility and VVIX satisfy a linear relationship, then the property of this stochastic volatil¬ 
ity follows from VVIX naturally. This relation can be seen as a benchmark to see whether 
the estimated stochastic volatility factor is accurate enough and also provide the empirical 
evidence for its model specification. From the historical data of VIX and VVIX, we find 
both of them are mean-reverting. Furthermore, through a formal test we find there exists 
evident co-jump between them. Thus we conduct an empirical analysis of the double-jump 
stochastic volatility model and its nested models using Markov-Chain-Monte-Carlo (MCMC) 
method with historical data of VIX and VVIX. We provide evidence consistent with jumps 
in both VIX and the volatility and that jump intensity is stochastic. We demonstrate the 
superiority of our main model through residual analysis and simulation result. 

The structure of this paper is as follows: Section reviews the related literature. Section 
shows the linear relationship between VVIX index and the stochastic volatility. Section 
analyzes model specification and sets ip our model. Section gives our empirical method. 
Section 1^ describes the data of VIX and VVIX that we use in this paper. Section [^summa¬ 
rizes the estimation results and provide empirical analysis. Section [^ concludes. 

2 Literature review 

To build a parameterized stochastic model for VIX, there are usually two different starting 
points. One is first modelling a multi-factor stochastic volatility process for S&P500 index 
and then derive a calculating formula for VIX index under this circumstance. For more 
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about multi-factor model setup, see e.g. Duffie et al. (2000), Gatheral (2008), Egloff et al. 


(2010), Gout and Kokholm (2013) and Papanicolaou and Sircar (2014). From the model 


assumption, the final VIX may be a combination of one or more factors (see e.g. Ait-Sahalia 


et al. (2014), Song and Xiu (2012), Luo and Zhang (2012), Lin and Ghang (2009)). The 


other method is to model VIX index directly which often has mean-reverting property. In 
this case, there are usually two ways to deal with the dynamics of VIX: affine and non-affine 


(see e.g. Mencia and Sentana (2013), Kaeck and Alexander (2013), Goard and Mazur 


(2013)). In the catalogue of non-a ffin e models, modelling logarithm of VIX directly is most 
popular and has been proved empirically better than affine assumption of VIX. 

Recently, modelling VIX with an additional stochastic volatility calls more attention. 
Along the literature, the existence of stochastic volatility in the dynamics of S&P500 has been 
widely proved and accepted empirically. Similarly, in terms of VIX index, the justihcation 


for stochastic volatility is also tested and verihed (see e.g. Wang and Daigler (2012), Huang 


and Shaliastovich (2014|)). Various authors has built VIX model with stochastic volatility 


and found the model with stochastic volatility factor outperforms the others without it. 


In Mencia and Sentana (2013), they make comparison of different models for VIX with 
VIX, VIX futures and VIX options as the data source. They use extended Kalman filter to 
estimate stochastic volatility and conclude that modelling logarithm of VIX with stochastic 
mean and stochastic volatility (named ’GTOUSV’ in their paper) is the best among all of the 
candidate models. In that paper, they hrst put forward the stochastic volatility of volatility 
model of VIX and model the volatility factor using pure jump OU process. [Kaeck and 


Alexander (2013) employ VIX data of nearly 20 years to estimate the VIX models with 


and without stochastic volatility using MGMG method. They model the volatility factor as 
a square-root diffusion process and prove the stochastic volatility model fits better for the 


VIX historical data. Barndorff-Nielsen and Veraart (2013) derive probabilistic properties of 


a class of stochastic volatility of volatility models. 
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From Huang and Shaliastovich (2014), we conclude that using only VIX index to infer 


the dynamics of the stochastic volatility factor of VIX is not a good idea. At least some 
additional data set must be taken into consideration. In terms of estimation, the data source 
matters. Estimation using various data sources can produce distinct empirical results which 


show their different information content (see e.g. Bardgett et ah (2013), Chung et ah 


(2011)). In the framework of parameterized SDE model for VIX, employing VIX and its 
derivatives (futures and options) has been implemented before. However the formulas for 
VIX options are generally vey complex and usually involve inverse Fourier transform which 
may increase the calculation burden. Closed-form solutions are thus not easy to obtain for 


VIX options. In that case, expansions is a good way (see e.g. Li (2013), Xiu (2014)). The 


information from VIX options are in essence equivalent to their implied volatility as other 
variables like time to maturity and strikes are known. As a volatility index for VIX options, 
VVIX just serves as the role of the implied volatility of VIX. Using VVIX index to estimate 
dynamics of VIX is still a vacuum in the literature and we will hll this gap in our paper. 

3 VVIX as a proxy for volatility of VIX 

3.1 A motivation 


In Kaeck and Alexander (2013), they put forward a stochastic volatility of volatility model 


for VIX. They assume the logarithm of VIX has a normal jump while its stochastic volatility 
factor satishes a square-root diffusion model. Denote by Y (t) the VIX index, we write this 
model under P, 


dY{t) = Kv{0-Y {t))dt+(t) + J^dN (t) 
du{t) = - u (t)) dt + a^\/u {t)dW^ {t) 


( 1 ) 


where i^dWy (t) ,dW^ (t)) = pdt. N (t) is Poisson process with constant jump intensity A, 

J? ~ {i‘", P,Y) 
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They use VIX data of 20 years as the input to estimate ([^ and give the estimated latent 
stochastic volatility u (t). A question arises naturally: how to guarantee and measure the 


accuracy of the estimation of this unobserved spot volatility? Bardgett et ah (2013) hnd 


that both of S&P500 options and VIX options contain some information that are unspanned 
by their counterpart. VIX is a volatility index summarized from S&P500 options while its 
diffusion part must reflect the implied volatility of VIX options. So employing VIX index 
as the only input to infer the dynamics of its stochastic volatility is doubted. While VIX 
options can provide wanted extra information, the computation burden rises quickly and 
sometimes complicated inverse Fourier transform is needed out of the complex expression 
of the price of VIX options. However, we will show in the next part that there exists a 
simple linear relationship between the stochastic volatility factor and the VVIX index which 
is compiled from a strip of VIX options mentioned above. 

3.2 Linear relationship between VVIX and volatility of VIX 


Following Mencia and Sentana (2013) and Kaeck and Alexander (2013), in this part we 


build model on the logarithm of VIX instead of VIX directly. We show that if logVIX mean 
reverts to a constant central tendency with stochastic volatility and jumps in logVIX and 
volatility, then there exists a linear relationship between VVIX index and this stochastic 
volatility factor of the logVIX. This relationship can provide a gauge to see whether a 
stochastic volatility model for VIX can be reliable. The similar idea of hnding a proxy for 


some unobservable factor can also be found in Ait-Sahalia and Kimmel (2007), Duan and 


Yeh (2010| and Ait-Sahalia et al. (2014). In this paper, VIX index is the underlying 


asset so its dynamics are observed under P measure. As VVIX is compiled from VIX options 
which is calculated under the pricing measure Q, all of the derivation involving VVIX below 
will be implemented under Q measure. 

Let Y (t) = logVIX (t) and assume that under Q, Y (t) and oj (t) follow a general affine 
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jump diffusion model 


dY{t) = kv{ 9-Y {t))dt+^/Uii)dW^ {t) + J^dN (t) 

du (t) = (a^ - K^oo (t)) dt + aaj^/ooJt)dW^ (t) + J^dN (t) (2) 

where we assume (^dWy (t),dW^ (t)^ = pdt. N (t) is a Poisson process with stochastic 
jump intensity A {t) = Aq + Aica {t) at time t for analytical treatability. The jump for VIX 
and its volatility factor are characterized by 


j: 


Q 


Y 




N K «) 


J3 ' N(P,(aif 


Similar to the idea that regarding VIX square as the expectation of quadratic variation 
of the logarithm of the S&P500 index under the pricing measure approximately (see, e.g. 


Ait-Sahalia et al. (2014)), we set 


VVIXly, = - 


Ef> (^jX(s)d,^+E3\^'£AYHB) 

With simple calculation from (|^ we obtain 


( 3 ) 


e: 


Q 


r-i+T 


u (s) ds I = 




Q 




U! (t) + 


T 




Q 

EX 






Ai/ici 


— oiQUj (t) + (3 q 


( 4 ) 


and 



+ i^y) ) i^o + (' 5 )) 

+ i^y) ) ('^o'T + EfXq + XiOqU (t)) 


( 5 ) 
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Combine ([^, (|^ and we finally have 

1 r 


VVIXl^, 


T 

1 

r L 


agu (t) + /3q + ( ifJ-y)^ + {(Jy) ] {XqT + Ai/?q + Xiagu (t)) 


I^Q + (+ i^y) ) (Ao't + XiPq) ) + f 1 + Ai f {fXyY + {o"y) ) ] aquj (t) 




A{t) + B (t) u (t) 


( 6 ) 


Relationship (|^ can be seen as a benchmark for the estimated volatility factor. The 
dynamics of VVIX can reflect the property of uj {t) more directly than the indirect impact 
of VIX option. It can provide more intuitive empirical evidence for the model specification 
for u (t) which will be seen in Section 

3.3 Examination using the benchmark 

If we define a suitable set of risk premia specification for this model to guarantee that its Q- 
counterpart remains the same structure as that under P, then the linear relation in Section 


3.2 holds. We want to use VVIX index as a benchmark (let = 0, Ai = 0) to judge whether 
the estimation for u {t) from only VIX index is reliable and reflect the real evolution of the 
market. As the VVIX data only starts from 2007, so in this part, we use VIX data from 


Jan 2007 to Sep 2014 to estimate this model again using the same method as Kaeck and 


Alexander (2013). We plot the estimated volatility factor u (t) and VVIX time series of the 


same period in Figure The correlation between this posterior volatility and VVIX index 
is only 0.4193. Although some of the peaks of estimated u {t) coincide with VVIX, more 
inconsistence between them appears. This indicates that when we use VIX index as the 
only data source to sample latent variable oj{t), it could only provide limited information 
about the dynamics of its stochastic volatility and the stochastic volatility is unspanned by 
the VIX in some sense. To obtain more accurate uj (t), the relationship between it and the 
VVIX index can be utilized. 
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4 Model specification and setup 

4.1 Model specification 

The log-normal Ornstein-Uhlenbeck model is put forward by ?. Ever since, modelling the 


logarithm of VIX or VIX futures is considered in Psychoyios et al. (2010) and Huskaj and 


Nossman (2013). Mencia and Sentana (2013) and Kaeck and Alexander (2013) compare 


and examine different model specification for VIX dynamics. Both of them conclude that the 
setup for modelling logVIX as an affine jump process is superior to modelling VIX directly 
which is consistent among all of the model specifications in more detail. So in our model, 
we also study the affine property of log VIX. 

Since there exists such linear relationship between VVIX and u (f), the VVIX index can 
be seen as a proxy for this unobservable variable. The joint modelling of VIX index and its 
stochastic volatility is thus equivalent to the joint modelling for VIX index and VVIX index. 
In this sense, the model should reflect some of their joint property. 

Both VIX and VVIX have the mean-reverting property. For VIX, it may mean revert to a 


constant or stochastic central tendency. In Mencia and Sentana (2013), they make both as¬ 
sumptions and examine the model performance respectively. As their data source consists of 
VIX, VIX futures and VIX options, the stochastic central tendency of VIX models performs 
better. In fact the specification of the central tendency of VIX is mainly characterized by the 
information from VIX futures while the VIX options play a relatively light role. However, 


from the derivation in Section 3.2 we know that the expression of VVIX index is irrelevant 
of the drift part of VIX. In fact, this is consistent with its stochastic volatility role. As our 
paper mainly concerns about the impact of VVIX data for the estimation, we make a simple 
assumption that the VIX mean reverts to a constant central tendency. As VVIX has the 
similar empirical property, we also assume the stochastic volatility of VIX has a constant 
central tendency. 

From historical daily data of VIX and VVIX from Jan 2007 to Nov 2014, we observe 



















that there exists evident co-jump between the two index, no matter positive or negative 
jump happened. To make a formal test to verify this phenomenon, we adapt the method 


m 


Bollerslev et ah (2008) to the lower sampling frequency (see also Gilder (2009) for the 


practice of this method for daily data). The testing procedure is divided into two steps: 
hrst, we show there exist jumps in both VIX and VVIX and second, the VIX and VVIX 
have common jumps. 

To carry on the hrst step, we assume a process X (to be VIX or VVIX) is observed in 
[0, T] at daily times t = 0,1,..., T and denote the time series by Xt, t = 1, 2,..., T. The 
return process rt = Xt — Xt_i,t = 1,2,... ,T is also dehned. We compute the n-day rolling 
sample estimates of realized volatility. 




k=0 


t—k 


and bipower variation 


BV, = ^ 


E n ,—1 


The relative contribution measure 


RJf = 


RVt - BVt 


(7) 


( 8 ) 


RVt 


(9) 


follows from (Q and ([^ immediately. The tripower quarticity for daily changes is defined by 

( 10 ) 




where p. 4/3 = 2^/^r T Finally, the statistic 

RJt 


Zt = 


TPt 

RV? 


( 11 ) 


y [( 7 r/ 2 )^ + TT - 5] 7max(^l, 

is constructed using (0), (0 and ( |To| to test whether a jump occurs at day t. We reject the 
null hypothesis of no jumps at a% conhdence level il\zt\ > where <F is the cumulative 

normal distribution for a given day t. 
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To implement the second step, denote the VIX and VVIX by and Assume VIX 
and VVIX index are observed in [0,T] at daily times t = 1,2,... ,T, the time series are thus 
XI, t = 1,2,... ,T,i = 1,2 respectively. Given the return processes rj = XI — Xl_i,t = 
1,2,... ,T,i = 1,2, we calculate the contemporaneous correlation 


cpt = J2 


n—l 

k=0 


r 


1 

t—k 


r 


2 

t—k 


and study the studentized statistic 


^cp,t 


cpt - cp 


( 12 ) 


^cp 


where 


and 


"P=r-(n-l) 


^cp - 


V {cpt - cpf 

< Jf—'n. 


1/2 


\T — {n — 1) ^t=r> 

at time t. We reject the null hypothesis of no common jumps at a% conhdence level if 
\zcp,t\ > ‘^i"-a /2 where <h is the cumulative normal distribution for a given day t. 

Employing the methods given above, we test the jump behavior of VIX and VVIX from 
January 3, 2007 to November 26, 2014. Given the 5% signihcant level, 222 days for VIX and 
141 days for VVIX out of 1939 days indicate the signihcant jump for the hrst step. In the 
second step, 131 days call for co-jump. Thus the specihcation for co-jump in VIX and VVIX 
is justihed and this phenomenon provides an important foundation for our model setup. 


4.2 Basic model 


This demonstrates that in addition to the diffusion part, we should assume jump in both Y (t) 
and u (t) and the jumps should be dominated by a single Poisson process. The jump intensity 
may be constant or state-dependent on the affine factor. In this paper we assume it is affected 
by u {t). The assumption for constant or stochastic jump intensity will be examined below. 
As both positive and negative jumps appear, we make the normal distribution assumption 
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for the jump size. For logVIX this may be a sensible assumption. While for the square root 
diffusion plus a jump for cu {t), as jump is a rare event for the historical path, the assumption 


is also acceptable. For more previous work on jumps in volatility, we refer to Duffie et ah 


(2000) 

Eraker et al. 

(2003) 

Eraker 

(2004) 

Todorov and Tauchen 

(2011 

) and 

Amengual 


and Xiu (2014). 


We thus assume that under Q, 


dY{t) = Kv{e-Y {t))dt+^/^)dW^ it) + J^dN it) 
du (t) = [a^ - K^uj (t)) dt + aoj\/u (t)dW^ (t) + J^dN (t) 

where (^dWy (t ), dW^ (t)^ = pdt and N (t) is a Poisson process with stochastic jump inten¬ 
sity A (t) = Ao + Aicu (t) at time t. 


j: 


Q 


Y 




We specify the risks of price between Q and P about Brownian motions as 


dW^{t) = dW^ {t) 
dW^it) = dW^ (t) - 


then under P, 

dY (t) = [kv {6 — Y (t)) — (^)] dt + \/00 {t)dWY (t) + JydN (t) 

duo it) = {a^ - K^oj (t)) dt + aujY'oj (t)dW^ (t) + J^N (t) (13) 

where (dWy (t ), dW^ (t)) = pdt and N (t) is a Poisson process with stochastic jump inten¬ 
sity A (t) = Ao + Aicu (t) at time t. is the speed of mean reversion under P. 
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The jump sizes are characterized by 



The parameter set under P is denoted by 

which the parameter set under Q is summarized as 

©M 1 } 


We also assume that there exists a pricing error for VVIX from our theoretical model, 

s.t. 

VVIXl,^ = A(T) + B(T)oj(t)+E (14) 

and 

£ ~ (O, CTp) 


we need to estimate 


©£ — Wp} 


We call the general model (13) the SVJJ-S model (stochastic A). If we let Ai = 0, it 


reduces to the SVJJ-C model (constant A) model. If we further let (J^) = 0, it collapses 
to the SVJ-C model (constant A) model. Finally, when there are no jumps, i.e., Jy = Jui = 0, 
we call it SV model. We want to examine these models using the real market historical data 
of VIX and VVIX to see: 1, whether adding the jump into VIX and u {t) can improve the 
VIX model signihcantly; 2, whether the jump intensity is constant or stochastic. 
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5 Model inference with VIX and VVIX 


In this part, we use VIX and VVIX index data from January 3, 2007 to November 26, 
2014 to estimate the models. In total we have 1991 daily observations for VIX and VVIX 
index respectively. We adopt MCMC method as the estimation method. Compared with 
maximum-likelihood estimation (MLE), generalized method of moments (GMM) and some 
other methods, MCMC has two advantages that adapts to our aim. First, not only does 
MCMC estimate the unknown parameters, it can also provide posterior estimated latent 
variables such as stochastic volatility, jump times and jump sizes. These variables are fun¬ 
damental and essential for subsequent empirical analysis and model comparison. Second, 
MCMC is very efficient for implementation. For more details about applications of MCMC 


method in hnance, we refer to Johannes and Poison (2003), Eraker et ah (2003) and 


Amengual and Xiu (2012). 


Denote by the parameters set by © = (0p, ©m, ©e), the latent variables by Z and the 
observed data by Y = {VIX, VVIX), for some model M we are interested in the joint 
posterior of parameters and latent variables given data: 


p(©,Z|Y,M) ap(Y|©,Z,M) ■p(©,Z|M) 


We assume the market data are observed daily. Let the time interval A = 1/252 be 
one day, assume we have T observations Va , 0 < i < T -|- 1 for logarithm of VIX. A time 


discretization of the dynamics (13) with time interval A gives 


Va - Y(i_i)A = {kvO - Ey Y(i_i)A - ?va;(i_i)A) A -h -h j/a^^a 

(UjA ~ — ^i^l^(i-l)A) A J- (T^\J (Up-i)AAe^ J- (15) 


where and are correlated Normal variables with correlation p, and are normal 
with different parameters. 

Denote by Va = Va - JiA^iA ioi 2 < i <T + 1 and UiA = uJiA - iiA^iA for 2 < i < T, 
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then we transform from (15) to the jump-adjusted processes 


hiA — flO + fllh(i_i)A + 0-2^(i-l)A + \/ ^^(i-l)A'^efA 

= Co-|-CiCi;(j_i)A + Y^wprpA^e^ ( 16 ) 

where Oq = KyO^, Oi = 1 — 02 = —<tyA, cq = a^^A, Ci = 1 — k^A. In this part, we 

will apply YiA, 0<i<T-|-lto estimate latent variables 

i^iA, 1 < i <T 

^*A,jfA and 2 < i<T + l 

As the joint posterior distribution p (©, Z|M) are not known in closed-form, the MCMC 
algorithm samples these parameters and latent variables sequentially from posterior condi¬ 
tional distributions as follows: 

spot volatility: p (ul^2 \^Sa^ ^i~A\ jiA~^\ jiA~^\ Y j 

jump time : p l^iA > jiA~^^ > ^) 

jump size in VIX : p \n2l , , hA~^^ , , Y ) 

jump size in volatility : p \n22, i^Ia ) JiA^ > ^ Y j 

parameters : p ^0^^^ \n\2 , uil2 > JiA^ ^ JiA^ ^ 0(9-b, Y j 

where g represents the iteration times. In this paper, we sample 5000 times and discard the 
hrst 2000 samples. 

5.1 Estimation Strategy 

In this part we consider the sampling method for the latent variables and parameters. We 
will discuss the corresponding algorithms for the stochastic volatility Ut, the Q-parameters 
0M and the pricing error parameter 0^. For jump times, jump sizes and parameters in 0p, 
the sampling methods are standard and Appendix gives detailed algorithm. 
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Sampling the stochastic volatility Ut should take the information from both VIX and 


VVIX into consideration. Utilizing the linear relationship in (14), we use random-walk 


metropolis method to sample Ut. Let ''' >'^(f+pA’''' where 

the index [g) represents the iteration times, we specify the full conditional density as 


n-Ma-^) 0 ( 5 - 1 ) y 
Pi ^ ’ ^iA 5 JiA 5 JiA T ^ 


oc 


w; 


(5) 


exp 


lA 


■ exp 


exp 


{Cf^ + Df^-2pCi^Di^) 

2(l-p2) 


[vvixfi, - A(t) - B(T)4‘2y'' 


2a p 


^ii+l)A + -^p+l)A “ 2pC(i+l)A-D(i+l)A 
2 (l-p2) 


where 


and 


CiA = 


-DjA — 


yiA - jTa ^^^iA - do- aiY(^i-i)A - a2UJ^^^ 


(*-l)A 


a 


(i+l)A 


^(£i)A^ 
^!a - - Co - 


V ■5(5-l)„(5-l) n n V n , 4^) 

Lp+1)A 4(i+l)A’^(i+l)A airp_i)A 0,2UJiA 


di’A 


D 


(i+l)A 


, ,( 9 - 1 ) _ .■‘^(5-l)„(9-l) _ r — r ! 

^(i+l)A 2(i+i)A'yi+l)A ^0 ^l^iA 

,(5) 




for 2 < z < T — 1. The case for z = 1 and i = T follows similarly. Note that this target 
density contains information from both VIX and VVIX. 


The Q-parameters ©m are related to the observed VVIX index through (14). We thus 


use random walk metropolis method to sample these parameters with the target density as 


exp 


EL {VVIX^;a -A(t)-B (r) 

2al 


\/^ap 

For the pricing error parameter ©g = {ap}, conditional on VVIXf^ and UiA, CjA = 
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VVIXf^—A {t)—B (r) uJiA ~ N (0, dp). Assume the prior for Up is 71 ^ 2 ^ (up) ~ InvGam Q^cr |,2 

we then sample using InvGam {^*^2 ^■, 0 *^ 2 ^ with a*^ 2 ^ = ao-si + and 0*^22 = 

YA^2{vyiX^A-Xr)-B{r)u2i^f 
«<72 2 + -^-2-^• 

5.2 Model diagnostics and specification tests 

5.2.1 Residual analysis 


Given the sampled posterior latent variables (spot volatility, jump times and jump sizes) and 
parameters, we can construct several statistics to test and assess the ability of the model to 
£t historical data. Recall the discretization of Y [t) during the MCMC estimation 


YiA — flo + 0-lY(i-i)A + 0-2(^{i-l)A + a/ <^(i-l)A 


iA 


where Y^a = Y^a — for 2 < i < T + 1 and oq = nyOA, oi = 1 — nyA, 02 = —?vA. 

The representation of follows immediately and given by 


- 

^iA ~ 


YiA ~ Q-q — aiT(j-i)A ~ Q'2UJ(i-i)A 


,2<i<T + l 


(17) 


a/ A;(j_i)AA 

With the estimated variables and parameters at hand we can calculated these residuals 
immediately. We will compare the Q-Q plot of the residuals of different models. If these 
residuals follow standard normal distribution approximately, then model performs well for 
fitting historical VIX index. If there exists big discrepancy between the residuals and standard 
normal distribution, the corresponding model must have potential for further improvement. 
With 

<^*A = Co + cia;(i_i)A + a/ cap-i) a Ae^ 

where cUjA = i^iA — JiA^jA for 2 < i < T, and Cq = a^A, ci = 1 — k^A. We can calculate 
the residual for u {t) similarly 

l^iA — Cq — Cia;(j_i)A 


^iA — 




-,2<i<T 


>^/^{i-l)AA 

This residual can be used to compare the various model for u (t). 


(18) 
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The jump times of Y (t) and u (t) can be used to test whether the jump intensity is 
constant or stochastic. If the posterior sampled jnmp times are clustered, the constant jnmp 
intensity assumption is rejected. 

5.2.2 p—value method 

We also perform simulation study using the posterior parameters to test different specihca- 
tions. We hrst specify some statistics that can reflect the dynamics of VIX and calculate 
these statistics for the logVIX data. Then for every model, we simnlate many trajectories for 
Y with the same sample size as the VIX data using the estimated parameters from MCMC 
results. With the simulated trajectory, we calculate the sample statistics and compare them 
with that obtained from original VIX data. More specihcally, we use the following reference 
statistics 

• standard deviation 

• skewness 

• kurtosis 

• maximnm 

• minimum 

• maxjump: the highest positive changes in the index 

• minjump: the highest negative changes in the index 

• avgmaxlO: the average over the 10 largest positive changes 

• avgminlO: the average over the 10 largest negative changes 

• various percentiles of daily changes. The percentiles are denoted by percNUM where 
NUM indicates the percentage level. 
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Denote these statistics calculated from logVIX data by 0^, k = 1, 2, • ■ ■ ,10. Then, for 
a given model, simulate N trajectories for Y using the estimated parameters from MCMC 
results. For the nth simulated trajectory D, 1 < n < iV, calculate the statistics above which 
is denoted by fc = 1, 2, • • ■ ,10. For every /c, 1 < /c < 10, calculate 


Pk = 




N 


(19) 


where 1^ is the indicator function. Too high or too low P/c, 1 < A; < 10 indicates that the 
given model may distort from the genuine form. For more details about this method, we 


refer to ?? and Kaeck and Alexander (2013) 


6 Data 


In 1993, CBOE introduced VIX index and it serves as a benchmark for the volatility of the 
market. In September 22, 2003, the CBOE revised the calculation method of VIX which 
utilized a wider range of S&PSOO options and back-calculated the new VIX to 1990. The 
well-known generalized formula now for calculating VIX is 


VIX^ ((, = (A') 


Ft 

— - 1 
Ko 


T-t^ Kf T-t 

i ^ 

which utilizes a strip of OTM S&P500 options prices Q (Ki) and Ft is the forward S&P500 
index level derived from S&PbOO options. 

In March 14, 2012, CBOE released a new volatility of volatility index called VVIX. VVIX 
is a measure of volatility of volatility which represents the expected volatility of the 30-day 
forward price of the CBOE volatility index. The calculation method of VVIX is similar to 
VIX, it is calculated from the price of a strip of at- and out- of the money VIX options, i.e. 


VVIX^ ((, T) = (A') 


Ft 

— - 1 
Ko 


T-t^ Kf '' ” T -1 

i ^ 

where O (iPj) is the midpoint of the bid-ask spread for VIX option with strike Kt and Ft 
is the forward VIX index level derived from VIX option prices. Kq is the hrst strike below 
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the forward index level Ft. Using this method, CBOE has also calculated the VVIX index 
before the release data up to the start of 2007. We plot the historical time series of VIX and 
VVIX from Jan 2007 to Nov 2014 in Figure [Tj From the picture we can see that the range 
of VVIX is at a signihcantly higher level than that of the VIX. Like VIX, VVIX also mean 
reverts to its historical mean value which is nearly 80. Furthermore, they share some of their 
peak values, especially during the 2008 hnancial crisis. Compared with VIX, VVIX is more 
volatile and when VIX is high, the range of variation of VVIX widens. The statistics of VIX 
and VVIX from Jan 2007 to Nov 2014 are summarized in Tabled] 


7 Empirical results 


In this section we discuss the estimation result for VIX dynamics among different models. 
The parameter estimation for four models are summarized in Table and the simulation 
results are showed in Table For all of the candidate models, the estimates of p are positive 


and around 0.52 which is close to the result in Kaeck and Alexander (2013) (p = 0.659 
for SVJ model in their paper). As the parameter enters into the drift of the VIX, the 
estimation for 9 is relatively low compared to the mean value of log VIX market data during 
the same period, is signihcantly larger than kv and this rehects the fact that the volatility 
of VIX or VVIX is more volatile than VIX itself. 

Figure gives the estimated volatility process of VIX for four models. As we use VVIX 
index as a proxy for the volatility, this four processes present similar forms. The correlation 
between the estimated spot volatility and VVIX of the four models are 0.9781, 0.9782, 0.9822 
and 0.9807 respectively. The average posterior volatility for SV and SVJ models is slightly 
higher compared to the other two models. This can be explained that the addition of jumps 
in volatility reduces the demand on the volatility process. 

Figure]^ shows the Q-Q plot of the residuals of VIX calculated by (17) for all of the four 
models and Figure plot the time series form. From the upper left panel in Figure we 
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find the SV model is misspecified for it requires very large shocks to Brownian motion. This 
can also be seen from the upper left panel in Figure Compared to the other three models, 
the range of the residuals for SV model is signihcantly larger and there exists many large 
innovations. 

From the two upper panels in Figure we can see the tail of the residuals become 
slightly thin so SVJ model improves SV model better. Much of the big Brownian shocks can 
be absorbed into the jump part. The estimated jump size in VIX is reported in the upper 
left panel in Figure However, from the simulation results in Table we hnd that for both 
SV and SVJ model, there are one or more statistics whose p-value are out of the [0.05, 0.95] 


bound. In contrast, in Kaeck and Alexander (2013), they also test the SV and SVJ model 
(with normal jump) and demonstrate that all of the p-value are within the [0.05, 0.95] bound. 
This shows that the addition of VVIX as the poxy for volatility of VIX help detect the further 
space of improvement for the stochastic volatility of volatility model for VIX. We also turn 
to Figure [^calculated using (18) which compares the residuals of volatility processes among 
four models. The residuals of SV and SVJ models are evidently larger than SVJJ-C and 
SVJJ-S model. 

Next we come to the SVJJ-C and SVJJ-S models. As mentioned above, the residuals of 
volatility processes for this two models performs better than SV and SVJ models. This shows 
the impact of jump on the volatility. Figure [^describe the estimated daily jump probability 
for SVJ, SVJJ-C and SVJJ-S models. Evidently, with the jump in volatility added, the jump 
occurs a bit more frequently. We recall that for SVJ model in which no jump happens in 
volatility, the jump time are determined mainly by the information of the VIX index. While 
for SVJJ-C and SVJJ-S models, we sample the jump time using both information or signal 
from VIX and volatility (VVIX). As we assume that the jumps of VIX and its volatility 
factor are determined by the same Poisson process, a big jump in volatility may raise the 
jump probability. This means that not only does the volatility jump, but furthermore it 
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jumps more heavily than VIX. Return to Figure]^ the bottom panels for SVJJ-C and SVJJ- 
S performs better than the upper ones, this show that the jump in volatility can also have 
impact on the dynamics of VIX. The influence channel can be through moments of high 
order or extreme values which can be seen from Table [H 

Unlike transient Brownian motion shocks, the influence of jump in volatility is more 
persistent. After positive or negative jump, the volatility enters a new regime. As the 
diffusion part of VIX, its effect will last for a period. A simple empirical method for judging 
the existence of jump in volatility in some day is to compare the fluctuation of a period of 
VIX data before and after that day. For example, on Feb 27, 2007, the VIX jumped from 
11.15 to 18.31. Before this day for a long period, the VIX looked very tranquil with very 
negligible variation and stay around 11. However, after this turning point, the VIX became 
more volatile and large up or down occurred more frequently, ranging from 12.19 to 19.63 
during the next 20 days. In fact, on Feb 27 the VVIX index also jumped from 70.33 to 
110.42. If we calculate the average of VVIX index for 20 days before and after this day, the 
results are 72.54 and 96.55 respectively. This indicate that the volatility had changed from 
the original state to a new and higher regime and thus made VIX index more active. This 
effect cannot be achieved by just a single Brownian shock on volatility and should be caused 
by jump. Thus the SVJ model is misspecihed also from empirical observation. 

Figure [^describes the jump of volatility for SVJJ-C and SVJJ-S models respectively. The 
jump size are almost positive with only a big negative jump in SVJJ-S models. With the 
mean-reverting property, the volatility reduces to its mean level through negative Brownian 
innovations after a big positive jump. This also indicates that the impact of positive jump 
can be persistent and signihcant. 

From Figure]^ we observe that for SVJJ-C model, the jump times are clustered. This 
is extremely unlikely under the constant jump intensity assumption. We can also see from 
Table that the estimation of Ai in SVJJ-S model is signihcant above zero. These facts 
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indicate that the SVJJ-S model is superior to SVJJ-C model and depict the dynamics of 
VIX more accurately. When the stochastic volatility u {t) enters into a relatively high regime, 
more jumps happen and affect the dynamics of VIX. 

8 Conclusion 


This paper discusses the model specihcation for stochastic volatility models of VIX from 
information of VVIX. We construct a volatility proxy of VIX using VVIX index as the 
benchmark and study its role for improving the model assumption of VIX from empirical 
observations. Based on the joint behavior of VIX and VVIX we propose a double-jump 
stochastic volatility model for VIX. We use MCMC method to estimate and compare different 
nested models using daily data of VIX and VVIX. Based on this, we point out that the jumps 
in VIX and volatility are essential and statistically signihcant and analyze the impact of the 
jumps on VIX dynamics. We show the jump intensity is stochastic and state dependent. 
The use of VVIX brings the estimation of some Q-parameters. Compared with richer dataset 
composed of VIX futures and options, the accuracy of these parameters have potential for 
further improvement. The corresponding risk premia could be further specihed. This will 
be left for future work. 


Appendix: MCMC algorithm for inference 


The MCMC sampling methods for parameters under physical measure P, jump times and 


jump sizes provided here are standard. Our sampling algorithms here borrow from Johannes 


and Poison (2003), Kaeck and Alexander (2010) and Amengual and Xiu (2012). To set 


up, the jump-adjusted discretization of the processes under P are stated as follows: 


Va — Oo -|- aiY(i_i)A + Cl2^(i-l)A + \/l^(i-l)A 

= Co-f-Cia;p_i)A + 
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where Oq = KyO^, ai = 1 — KyA, 02 = —?yA, cq = ttojA, ci = 1 — k^A and I^a = 
hiA — jfA^iA, = A;iA — is to estimate P-parameters 

0P = {kv,(;v, 0, K^, ,af^, p,a^} 

and latent variables n^A, jTa JiAi 2 < i < T + 1. 

A. Sampling latent variables 

• Sample Jump times 


For niA,* = 2,3, • ■ ■ ,T, 


p(n,A = l|X, 0p,F) 

P ^JJa, A’iA |F(j_i)A, A'(i—1)A, "^iA J(i—1)A’J(i—1)A’ ' P (jT'iA 11 A’(j—1)A, F(i—1 )a) 

P ^h^A, A^jAlFi^j—1)A, 1)A, "^iA 1)A’J(i—1)A’ ' P {p'iA 'S|a'P_i)a, F(i—1 )a) 

where p (^F^a, a;iA|F'(j_i)A, cap.pA, n^A = -s, j(Vi)A’Jp-i)A> is a bivariate normal dis¬ 
tribution with mean 


Oo + aiF(i-l)A + fl2l^(i-l)A + S ■ 
Co + Cia;(j_i)A + s ■ 


and covariance matrix 


i^p-l)A 


A 


1 PC^U, 

.2 


p(^u, (y: 

and p (niA = Ijccp-pA, ^P-i)a) = (Aq + Aiccp.pA) A, 


For n(T+i)A, 


p(«(r+i)A = 1|A, 0p,F) 

P {X(T+i)A\YTA,^TATn{T+i)A = 1 Ata’JtA) ®p) ' P (^(r+i)A = IIAta) 
Ep(A"(r+l)A|hTA,C^rA,n(T+l)A = s, JpA^iTA^Qp) ' P (^(r+i)A = sEta) 

s=0 

where p (F(t+i)a|Fta, i^ta,''^( r+i)A = s-, Jtai JrAi^p) is a univariate normal distribu¬ 
tion with mean ao+aiFrA+a2i^TA+-s-jp’+i)A '^t’A A and p (n(T+i)A = IIFpa) = 
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(Ao + AiCUta) We can thus sample n^A for i = 2, 3, • ■ ■ , T + 1. 


Sample Jump Sizes 

When UjA = 1, we then sample using -A^ (x, where 

1 1 

A = 


UiA — Co — cia;(i_i)A /i, 
-D = -^^-h — 


JP 


and sample using A where 

1 1 


W, 




Al = 


B = 


^^p-l)A (1 - p2) A 

h^A — Clo ~ fllh(j_i)A — a2l^p-l)A ~ ^ (l^iA ~ Cq — Cia;(j_i)A — J^a) h 


+ 


^^p-i)A (1 - p2) A ■ 

When n^A = 0, the posterior distribution of j'^a JiA same with the prior 
distribution, i.e., 


P' -W (p'’. {‘’if) 


where the parameters are updated simultaneously. 
We thus sample jfA JiA ior i = 2,3, ■ ■ ■ ,T + 1. 

B. Sampling parameters 0p 

The parameter set under P is denoted by 


0P = {Kv,(^v,0,K^,fiy^,ni^,cT^,p,auj} 


• Sampling 6 

Assume the prior for 6* is A {pe, aj), we sample the posterior using 9 N{B/A, l/A) 


24 



with 


A 

B 


KyA 

2-^ nAT 


^9 ^ (1 -P^)^h-1)A 


he , 

^ a . ^ 


i=2 


YiA — Y{i-i)A + ^yh(i-i)AA + ^ya;(i_i)AA — pDi/\ ^ a;(i_i) a A 

(1 - p^)uJ{i-X)A 


• Sampling Ky 


Assume the prior for ny is N '^Ky) > sample the posterior using Ky rsj N{B/A, l/A) 
with 


A 

B 


a 


(0-yp-i)A)'A 

2^ n - 




pKv 


a. 


ky 


i=2 

T 

i=2 


(1 - p2)a;(i_i)A 

^iA — h(j-i)A + ?yi^(i-i)AA — pAa a/i^p- i)aA) (0 - F(,_i)a) 

(1 - p2)a;(i_i)A 


• Sampling <^y 


Assume the prior for ^y is N sample the posterior using ^y rsj A(5/A, 1/A) 

with 

^(i-l)AA 
1 — p2 

(^^iA — y{i-l)A — Ky {9 — A-Pa) a — pAaa/i^P-PaA^ 

1 -p2 

• Sampling k/) 

Assume the prior for is N j, we sample the posterior using N{B/AMA) 

with 


B = 


pKv 


a, 


-E 


w 


i=2 



A 

B 





(1 - P^) 


CJjA ^(i—1)A CIqjA CT^pCi/\ -sj A A 

(1 - p2) a2 


• Sampling (a^)^ 
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Assume the prior for these parameters are: ~ N (af jp ), ~ N ( afjp 

y \ l^y J \ t^oj 

~ InvGam 

Then we sample the posterior using 
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Table 1: Summary Statistics 


This table provides summary statistics for VIX and VVIX index from January 3, 2007, to November 26, 2014. 



Mean 

Volatility 

Skewness 

Kurtosis 

Min 

Max 

VIX 

21.9101 

10.3966 

2.1241 

5.7181 

9.89 

80.86 

VVIX 

85.9204 

12.8226 

0.8289 

1.0079 

59.74 

145.12 


Table 2: VIX Parameter Estimates 


This table shows the parameter estimation results for the four models using VIX and VVIX index data from January 3, 2007 
to November 26, 2014. Four each parameter, we give the mean and the standard deviation of the posterior. ”SV” denotes 
diffusion model with no jumps. ”SVJ” introduces jumps in VIX in the SV model with constant jump intensity, ’’SVJJ-C” adds 
double jumps in VIX and its volatility in the SV model with constant jump intensity. ’’SVJJ-S” assumes the jump intensity to 
be stochastic in the SVJJ-C model. 


SV SVJ-C SVJJ-C SVJJ-S 



Mean 

Stddev 

Mean 

Stddev 

Mean 

Stddev 

Mean 

Stddev 

Ky 

1.6800 

0.5733 

1.5765 

0.5600 

1.8611 

0.5688 

2.1093 

0.5866 

w 

-1.1869 

0.7718 

-0.8046 

0.7673 

-0.2702 

0.8305 

-0.1538 

0.7820 

e 

2.3500 

0.4404 

2.3090 

0.4446 

2.2704 

0.3954 

2.3312 

0.3120 


4.5162 

1.0284 

4.4308 

0.9973 

6.1132 

1.0650 

6.2849 

1.0645 


7.5104 

0.4314 

7.6866 

0.4676 

2.5996 

0.2584 

2.5674 

0.1958 


3.8549 

0.7807 

3.7683 

0.7540 

4.0781 

0.7882 

3.7938 

0.7308 

P 

0.5392 

0.0161 

0.5596 

0.0141 

0.5204 

0.0169 

0.4998 

0.0190 

CTu 

0.8560 

0.0724 

0.8207 

0.0117 

0.8848 

0.0853 

0.8461 

0.0372 




0.6550 

0.0492 

2.4295 

0.1787 

2.7557 

0.1332 

Ai 







1.6086 

0.1262 

Py 



0.1593 

0.0220 

0.1999 

0.0279 

0.1551 

0.0171 

Py 



-0.0520 

0.0037 

-0.0556 

0.0746 

-0.0960 

0.0306 




0.1075 

0.0172 

0.1121 

0.0132 

0.1231 

0.0108 

pi^ 





0.1872 

0.0226 

0.1430 

0.0239 






-2.0084 

0.0882 

-1.2046 

0.0547 






0.1307 

0.0165 

0.1420 

0.0161 

ap 

0.0599 

0.0077 

0.0592 

0.0071 

0.0563 

0.0082 

0.0612 

0.0076 
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Table 3: Simulation results 


This table reports the p-values calculated by for all the statistics of simulation results of VIX for different models. It 
describes the average comparisons of the statistics of historical data and the simulation paths from every given model. Very 
high or low p-values indicate the model’s inability to capture the VIX dynamics. ”SV” denotes diffusion model with no jumps. 
”SVJ” introduces jumps in VIX in the SV model with constant jump intensity, ’’SVJJ-C” adds double jumps in VIX and its 
volatility in the SV model with constant jump intensity. ’’SVJJ-S” assumes the jump intensity to be stochastic in the SVJJ-C 
model. 



Data 

SV 

SVJ-C 

SVJJ-C 

SVJJ-S 

stadev 

0.4493 

0.0670 

0.3351 

0.8593 

0.3733 

skewness 

0.9005 

0.1462 

0.0359 

0.2194 

0.7479 

kurtosis 

0.7806 

0.2453 

0.0150 

0.6250 

0.7186 

maximum 

4.7558 

0.0079 

0.1383 

0.6875 

0.4744 

minimum 

2.3984 

0.2418 

0.6040 

0.0769 

0.6896 

maxjump 

0.2267 

0.4897 

0.1016 

0.2805 

0.6358 

minjump 

-0.2422 

0.0953 

0.9539 

0.7698 

0.7509 

avgmaxlO 

0.1852 

0.5697 

0.1508 

0.4308 

0.2614 

avgminlO 

-0.1671 

0.3262 

0.9069 

0.7996 

0.5826 

percO.Ol 

-0.1345 

0.4271 

0.8488 

0.5005 

0.6311 

perc0.05 

-0.0931 

0.4504 

0.4957 

0.5048 

0.6277 

percO.95 

0.0928 

0.5604 

0.6182 

0.8245 

0.4581 

percO.99 

0.1370 

0.4478 

0.0264 

0.9379 

0.1622 
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2008 2010 2012 2014 


This figure shows the time series of VIX and VVIX index from January 3, 2007 to November 26, 2014. Both of them are 
mean-reverting and VVIX is at a significant higher level than VIX in terms of the range of values. 

Figure 1: VIX and VVIX index 
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The spot volatility in this figure is the estimated posterior volatility of logVIX in SVJ model with only VIX index as the data 
source. It shows the comparison of this volatility and the contemporaneous VVIX index 

Figure 2: Spot Volatility from VIX estimation vs VVIX 
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SV model 


SVJ model 




SVJJ-C model SVJJ-S model 




The figures show the estimated paths of posterior volatility oj (t) for four models. All of them are highly correlated with VVIX 
index. The level of the volatility in SV and SVJ models is slightly higher than that in SVJJ-C and SVJJ-S models. ”SV” denotes 
diffusion model with no jumps. ”SVJ” introduces jumps in VIX in the SV model with constant jump intensity, ’’SVJJ-C” adds 
double jumps in VIX and its volatility in the SV model with constant jump intensity. ’’SVJJ-S” assumes the jump intensity to 
be stochastic in the SVJJ-C model. 

Figure 3: Posterior volatility of VIX for each of the models 
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SV model 


SVJ model 



Theoretical Quantiles 



Theoretical Quantiles 


SVJJ-C model 


SVJJ-S model 




Theoretical Quantiles 


Theoretical Quantiles 


The figures show the Q-Q plot of the residuals calculated from each of the models using 0 with the estimated parameters as 
input. SVJJ-C and SVJJ-S models perform relatively better than SV and SVJ models. ”SV” denotes diffusion model with no 
jumps. ”SVJ” introduces jumps in VIX in the SV model with constant jump intensity, "SVJJ-C” adds double jumps in VIX 
and its volatility in the SV model with constant jump intensity. ’’SVJJ-S” assumes the jump intensity to be stochastic in the 
SVJJ-C model. 

Figure 4: Q-Q plot of the residuals 
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SV model 


SVJ model 



2008 2010 2012 2014 2008 2010 2012 2014 


The figures show the time series of standard innovations or residuals of VIX calculated from the estimated parameters using 
” SV” denotes diffusion model with no jumps. ”SVJ” introduces jumps in VIX in the SV model with constant jump 
intensity, ”SVJJ-C” adds double jumps in VIX and its volatility in the SV model with constant jump intensity. ”SVJJ-S” 
assumes the jump intensity to be stochastic in the SVJJ-C model. 

Figure 5: VIX Residuals 
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2008 2010 2012 2014 


SVJJ-C model 



2008 2010 2012 2014 


SVJJ.S model 



The figures show the time series of average jump sizes in VIX. ”SV” denotes diffusion model with no jumps. ”SVJ” introduces 
jumps in VIX in the SV model with constant jump intensity, ’’SVJJ-C” adds double jumps in VIX and its volatility in the SV 
model with constant jump intensity. ’’SVJJ-S” assumes the jump intensity to be stochastic in the SVJJ-C model. 

Figure 6: Posterior Mean of Jumps in VIX 
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SV model 


SVJ model 





2008 2010 2012 2014 


(D - 



2008 2010 2012 2014 



The figures show the time series of standard innovations or residuals of volatility of VIX calculated from the estimated parameters 
using . ”SV” denotes diffusion model with no jumps. ”SVJ” introduces jumps in VIX in the SV model with constant jump 
intensity, ’’SVJJ-C” adds double jumps in VIX and its volatility in the SV model with constant jump intensity. ’’SVJJ-S” 
assumes the jump intensity to be stochastic in the SVJJ-C model. 

Figure 7: Volatility Residuals 
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SVJ model 
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2008 2010 2012 2014 


SVJJ-C model 



SVJJ.S model 



The figures show the estimated jump probability of SVJ, SVJJ-C and SVJJ-S models. ”SVJ” introduces jumps in VIX with 
constant jump intensity and models the volatility using square root diffusion model, "SVJJ-C” introduces double jumps in VIX 
and its volatility with constant jump intensity. "SVJJ-S” assumes the jump intensity to be stochastic in the SVJJ-C model. 

Figure 8: Estimated jump times for SVJ, SVJJ-C and SVJJ-S models 
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The figures show the time series of average jump sizes in the volatility of VIX. "SVJJ-C” introduces double jumps in VIX and 
its volatility with constant jump intensity. "SVJJ-S” assumes the jump intensity to be stochastic in the SVJJ-C model. 

Figure 9: Posterior Mean of Jumps in volatility of VIX 
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